The Lame equation in parametric resonance after inflation 
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We show that the most general inflaton potential in Minkowski spacetime for which the differential 
equation for the Fourier modes of the matter fields reduces to Lame's equation is of the form 
V{<j)) = \(j} 4 /4 + K<f> 2 /2 + fi/(2(f) 2 ) + Vo. As an application, we study the preheating phase after 
inflation for the above potential with K — and arbitrary A, n > 0. For certain values of the 
coupling constant between the inflaton and the matter fields, we calculate the instability intervals 
and the characteristic exponents in closed form. 
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PACS numbers: 98.80.Cq, 98.80.Hw 



I. INTRODUCTION 



The phenomenon of parametric resonance plays a fun- 
damental role in the modern theories of preheating after 
inflation jl], ||, ||. At the end of the slow-roll phase, 
the inflaton field starts oscillating around the minimum 
of the potential in a coherent way. If this field is cou- 
pled to some other matter fields, the periodic evolu- 
tion can give rise to an explosive production of matter 
quanta thanks to the resonant amplification of vacuum 
fluctuations. This production is typically characterized 
by the exponential growing of the occupation number of 
those states whose momentum lies within certain reso- 
nance bands. This effect is essentially non-perturbative, 
and the corresponding production of particles is much 
more efficient than the traditional perturbative decay of 
the inflaton field during reheating. However, this non- 
perturbative character makes it extremely difficult to ob- 
tain exact results in most inflation models and, as a con- 
sequence, one must rely on numerical computations in 
order to obtain information about spectra or time evolu- 
tion of the particle production. Despite this fact, a few 
models are known for which it is possible to obtain the 
width of the resonance bands and the characteristic ex- 
ponents in an analytical fashion. In particular, the pure 
quadratic potential V{(j>) = m 2 </> 2 /2 for the inflaton field 
in Minkowski space yields the Mathieu equation for the 
Fourier modes of a scalar field \ coupled to the infla- 
ton as g 2 (f> 2 x 2 /2 (|, (§. In an expanding background, 
the quartic inflaton potential V(4>) = \<fi A /4 leads to 
the Lame equation for the corresponding matter fields 
modes This equation also appears in connection 

with some other combinations of these inflaton poten- 
tials, suchasF(0) = {<fi 2 — er 2 ) 2 ||. Parametric resonance 
driven by the latter potential has also been proposed as 
a mechanism for the generation of long-wavelength pion 
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modes from disoriented chiral condensates (DCC) 0. 

The relevance of the Lame equation in the study of 
parametric resonance stems from its unique analytical 
properties, which makes it possible to compute the reso- 
nance bands and the associated characteristic exponents 
in closed form ]|, ||. In this paper we extend the above 
results and derive the most general inflaton potential in 
Minkowski spacetime for which the corresponding matter 
modes equation reduces to the Lame form. The result- 
ing potential possesses a term of the form 4>~ 2 , appearing 
typically in the context of certain supersymmetric infla- 
tion models, see Refs. js| |). 

The paper is organized as follows. In Section II we 
present our main result, obtaining the most general po- 
tential for the inflaton field in a Minkowskian background 
leading to the Lame equation for the matter fields via the 
usual g 2 2 x 2 /2 coupling. For a certain choice of the pa- 
rameters of this potential, we explicitly determine the 
scaling factors necessary to reduce the equation for the 
Fourier modes of the matter fields to the Lame equa- 
tion. In Section III, we determine the range of values of 
these parameters compatible with the slow-roll approx- 
imation and the amplitude of the density perturbations 
observed by the Cosmic Background Explorer (COBE). 
For several values of the coupling constant g we compute 
(using the results of Refs. JE], ||) the resonance bands, the 
corresponding characteristic exponents, and the number 
density of particles produced. Our results show that the 
particle production in the preheating phase remains un- 
changed if a term fi/(2<fi 2 ) is added to the A</> 4 /4 model 
for a broad range of values of the parameter /i. We finally 
summarize in Section IV the main results of our paper. 



II. THE INFLATON POTENTIAL 

Let us consider the following Lagrangian density for 
the inflaton field 4> coupled to a scalar massless matter 
field x : 



C 



1 



V{cj>) 



3,2 2 

y0X ■ 



•odriguc@cucmos.sim.ucm.ef 



(1) 



2 



We will assume our background fields to be homoge- 
neous and isotropic, i.e., the space-time metric is of 
the Friedmann-Robertson-Walker form ds 2 = dt 2 — 
a 2 (t) dx 2 , and the inflaton field <f>(t) depends only on 
time. The corresponding classical equations of motion 
for the different fields are 

4> + 3H(j>+V'(<t>) =0 (2) 

X + 3Hx- ^d t d lX + g 2 <t?X = , (3) 
or 

where H — a/a, and the prime denotes differentiation 
with respect to the argument of V. Note that we are 
neglecting back reaction effects from the x fields in the 
inflaton equation (|^) assuming that the expectation value 
(x 2 ) is negligible during the first inflaton oscillations. 

Let us rescale the field as ij> = a 4>, x = a Xi an d work in 
conformal 77 time defined by a drj = dt. ft is then possible 
to rewrite the equations of motion as 



under a change of scale 

t = rx, <Kt) = fa f(x) , Xk(t)=X k (x). (10) 

In the latter equations m = k 2 € (0, 1) is the square of 
the modulus of the elliptic sine function sm = sn(cc, k). 
Comparing Eqs. (|J) and we see that we can take 
without loss of generality 



f(x)=±((3- 



2 \V2 

sn x) 



and 



£ = r 2 {k 2 + (3g 2 4>l), mn(n + l) =T 2 g 2 (j) 2 . 

Note that j3 > 1 for / to be real for all x £ R. 
It follows from Eq. M) that 



(11) 

(12) 
(13) 



d 2 4> _|_ ffl 3y'/ -i^ 1 < ^ a 
dr\ 2 a drj 2 



= 



77 - o l d l x 
dr) z 



2 ~ 2 A 1 d 2 a A 

9 <P X 77X = 0- 

a drj z 



(4) 
(5) 



In the particular case in which the potential is V((f>) — 
A0 4 /4, it can be shown [[| that the scale factor grows as 
0(77) oc 77, and hence the last terms in Eqs. (f|)-(||) vanish, 
while the term a 3 F'(a~V) reduces to A0 3 . The equa- 
tions of motion are thus expressed in Minkowskian form 
in terms of the rescaled fields. However, this is not true 
for an arbitrary potential, and for that reason in the fol- 
lowing we will concentrate only in the Minkowskian limit 
(a = 1). This can be considered as a first approxima- 
tion to the full problem and, in particular, it will allow 
us to obtain analytic results in some cases. Nevertheless, 
we expect that our results will carry over to an expand- 
ing universe with minimal quantitative changes (though 
rather important analytic differences) when the potential 
differs from A</> 4 /4 by a small perturbation. 

Let us then consider the equations of motion in 
Minkowski space-time 



fx 

dt 2 



didix + g 2 <j> 2 x = 



(6) 
(7) 



We shall now derive the most general inflaton potential 
V{4>) for which the differential equation for the Fourier 
modes of the matter fields 



d 2 Xk 
dt 2 



+ (k 2 + g 2 4?) X k = 



(8) 



reduces to the Lame equation 



dx 2 



(£- mn{n + 1) sn 2 x) X k =0 (9) 



where the integration constant E is the energy density of 
the inflaton field. Using the well-known identities 



dx 



sn x = cn x dn x 



and 



ca 2 x=l 



P + f 2 {x), dn x = 1 - 777/3 + to/ 2 (x). 

it is immediate to show that 
2 



Id.) 



(14) 



with U and Ef given by 
to/ 4 1 



2/2 



E f = -- [3to/3 2 - 2(1 + m)(3 + l] 



Comparison of Eqs. ( |l3| ) and (14) using ( |l0| ) shows that 
the inflaton potential V(cf>) must be of the form []l6| 



VW=\Z r + KZ- + £ s + V 



4 ' ~ 2 ' 2(f> 2 



(15) 



where the potential coefficients (A, K, /i) and the en- 
ergy density E are related to the parameters (/?, m,T,fa) 
through the following equations: 



A = 



2to 

4 



K = — (l + m-3m/3) 



i x = ^p{p-l){l-mP), 



E 



ii. 

2r 2 



[3to/3 2 -2(1 + to)/?+1] +V . 



(16) 

(17) 
(18) 



Note that the positivity of A forces n to be non-negative 
for V to have a global minimum. We have thus shown 



3 



that, in order for Eq. (ph to reduce under scaling to 
Lame's equation @, the infiaton potential V(4>) must 
necessarily be of the form (|l5|). 

One must now prove that any potential of the form 
( |l5| ) leads to Lame's equation (|9|). For this to be the 
case, Eqs. (|l^)~(|l8|) should have a solution for arbitrary 
values ofA>0,-ftT,/_i>0 and E > V^ lin where V m i n 
is the minimum value of the potential ( |l5[ ) (the value 
E = Knin must be discarded, since then cj) reduces to the 
trivial constant solution cf) — m i n ). We shall verify this 
only for the case K = 0, which we shall study in detail 
in the following sections. 

For K — the second Eq. (|l6|) can be used to solve for 
(3 in terms of to, obtaining 



The condition (3 > 1 then implies that 



< to < - . 
~ 2 

When < to < 1/2 the remaining equations (|T7|)-(|T8|) 
are easily solved, with the following result: 



r 2 = ±{\^y^D{m) 



E = Vi 



TO — TO + 1 

D 2 (m) 



/ix 1 / 3 3to 
x) D(m) 



(20) 
(21) 



where 



D(m) = 



(1 + to)(2 — to) ( — — to 



-l 1/3 



and Vi = l(Xfi 2 ) 1/3 is the minimum of the potential fll5| ) 
for K = V = 0. Since 



E(m) 



m — to + 1 
D 2 (m) 



(22) 



grows monotonically in the interval m 6 (0,1/2), with 
i?(0) = 1 and lim„ wl / 2 E(m) = +oo, it follows that 
Eq. ( |2l| ) uniquely determines to £ (0,1/2) for arbitrary 
values of E > V± + Vq = V m i n . Eq. ( p0| ) then yields r and 
cfo for arbitrary A,/i > 0. For to = 1/2 Eqs. ©-© 
simplify to 



A : 



M = 0, 



E = 



At 2 



These are the equations obtained in Ref. || for the A(/> 4 /4 
potential, which again uniquely determine r and in 
terms of A > and _E > V m i n = 0. In fact, for to = 1/2 
Eq. (|l|) implies that @ = 1, so that (0) yields / 2 (x) = 
cn 2 (x, l/\/2). 



III. EXACT RESULTS FOR THE A 4> 4 /4 + u/(2(f> 2 ) 
MODEL 

Potential terms of the form [i(j)~ p have been consid- 
ered in different contexts in the literature. They appear 
in hybrid inflation models jl(j, or in the so-called inter- 
mediate inflation jO]. But it is probably in the context 
of supersymmetric models of inflation where this kind of 
terms is more relevant. In fact, these contributions arise 
generically due to non-perturbative effects in supersym- 
metric gauge theories |T^], the scale /x and the exponent 
p depending on the particular gauge group in the theory. 
Inflationary models based on this kind of potentials have 
been studied in Refs. || 

Motivated by the connection with the Lame equation 
established in the previous section, we shall consider an 
effective potential during inflation and preheating of the 
form 



(23) 



We shall derive an analytic expression for the instabil- 
ity intervals and the characteristic exponents for certain 
values of the ratio g 2 /X. In principle, higher-order terms 
could appear in the full potential, but we will assume 
that their effect is negligible during inflation and pre- 
heating. We could also have considered the contribution 
of a mass term as shown in the previous section, but in 
order to obtain closed expressions we will ignore it in 
what follows. 

The minimum of the potential is placed at </> m in = 
(///A) 1 / 6 = a, while the constant Vq must be taken as 
Vq = —Vi, i.e. Knin = 0, for the cosmological constant to 
vanish. The infiaton field <j) oscillates around m i n with 
amplitude a U/l + m — y/T— 2to) /(2y/ D(m) ) and pe- 
riod 2tK(k), where K (k) is the complete elliptic integral 
of the first kind Note that, in contrast to the pure 
X(f> A /4 model, <fi m i n is nonzero and the oscillations are not 
symmetric about (f> m - m (see Fig. [l]) . 

Let us introduce the following notation: cf> — <p/a and 
a = y/na/Mp. The slow- roll parameters e and 77 for this 
model are given by: 



Ml 
16tt 




and 



M 2 P V" 



1 



2a 2 



30 4 + 2c)) 2 



(24) 



(25) 



The end of inflation occurs when the slow-roll approx- 
imation breaks down, i.e., for e ~ 1 or 77 ~ 1. In the 
case a <C 1 the value </> cn d at the end of inflation is 
given to a very good approximation by <^ on d — \/3/2 a -1 . 
However for a 3> 1 the corresponding value behaves as 
(f) cn d ~ 1 + 5 a -1 . The initial value 060 corresponding 



4 



<f>/M f 




1CT 6 M P t 



log 10 a 



FIG. 1: Oscillations after inflation of the inflaton field (j> (in 
units of Mp) as a function of time t (in units of W 6 /M p ) 
computed from the exact formula ([h)|)-([t~l|). We have taken 
A = 9.091 ■ 10~ 13 , fi 1/6 = 5.553 • KT^Mp, m = 0.417, corre- 
sponding to a = 1; see the discussion in the text for details. 
The horizontal axis is placed at (j>min/M p = 1/yfn. 



to 60 e-folds before inflation ends is determined by the 
condition 



60- - 



8tt 



^>60 



V 
V 



7 



(26) 



It follows that 060 — \/123/2a _1 if a <C 1 and </>60 — 
1 + 5 _1 if <5 3> 1. Accordingly, the amplitude of the 
density perturbations at 4> = 4>qq can be written as: 



FIG. 2: log-log plots of 10 12 A and /i 1/6 /M P versus a for 
10" 3 < & < 10 3 . 



FIG. 3: Value of the square modulus m as a function of a. 



S H (k) = f 



2 

We thus get S H {k) 



2tt V 3 ' 2 



3 M%\V'\ 



(27) 



-s(* 
a 6 — 



30| 



3/2 



/41A 



' 123 

v 5 



4o 



96 ~4 



41 



a 



) for a <C 1 and 
2A<5 for 5 ^> 1. The COBE normalization 



Sji(k) — 5 • 10 5 determines A (and thus ji) as a function 
of a, see Fig. ^. In the limit a = we get A ~ 10" 12 , 
while A ~ 5 • 10" 12 6r 2 for <5 > 1. Since /_t 1/6 < M P we 
get A > 2 • 10~ 18 when a ^S> 1. This in turn implies that 
3 < 2 • 10 3 . 

The value of m is obtained from Eq. (^l|), with the 
energy density of the inflaton field given by 



E = V(0 ond ) = V t 



30 



2 

end 



3^ 2 „d 



(28) 



For a <C 1, to a very good approximation we have E 



3 ^, 



a . From Eq. (E2J) we thus get m 
the other hand, for d> 1 we have E ~ 



Jo" 



On 

<5~ 2 , leading 



In Fig. H we present the plot of m as a 
function of a for a between and 10. 



to to ~ | (5 



If 77 is a positive integer, the Lame equation (^) pos- 
sesses exactly ti + 1 instability zones as the parameter 
£ takes values on the real line, whose corresponding so- 
lutions grow exponentially at either ±oo |l4[ . However, 
only those instability bands for which the squared mo- 
mentum k 2 given by Eq. ( |l2| ) is non-negative are phys- 
ically significant. According to Floquet's Theorem, the 
solutions of the Lame equation in an instability band can 
be written as Xk(x) = e ±AIfe2: Pk{x), where Pk{x) is a pe- 
riodic function and the characteristic exponent [ik has a 
nonzero real part. The occupation number for particles 
with momentum k produced in the preheating phase can 
be estimated in terms of the characteristic exponent as 
Nk{t) — e 2AIfct / r IJ. Remarkably, for integer n the char- 
acteristic exponent for the Lame equation is given by an 
exact formula involving a quadrature which, at least for 
the lowest values of n, can be expressed in terms of ellip- 
tic integrals ||, ||. We shall omit here most details and 
quote from the above references only the main steps for 



5 



the derivation of this formula. 

The key ingredient for obtaining the formula is the 
construction of an exact expression for the product of 
two linearly independent solutions of the Lame equation 
jl7j . The product of two such solutions satisfies the third- 
order differential equation given by 

2p{z)M"'{z) + 3 P '(z)M"(z) + \p"(z) + 2(8 

+ mn(n + l)(z-l)) M'( z) + mn(n + l)M(z) = , 

(29) 

where z = cn 2 x, and 

p(z) — (1 — m)z + (2m — l)z 2 — mz 3 . 



If n is a nonnegative integer, Eq. (|29j) is satisfied by a suit- 
able polynomial of degree n, which (following ||) shall 
be written as 



M (n) (z) 



i=0 



(n) 

with the normalization condition <z = 1. Since z is a 
periodic function of x, the polynomial Mt n \(z) coincides 
with the product of two linearly independent solutions of 
the Lame equation (^) when £ belongs to an instability 
zone. It was shown in Refs. || |(| that the characteristic 
exponent is determined by the definite integral 



(30) 



C, 



(n) 



2K( K )J JpJz) M{n) {z) 



(31) 



where 

C 2 n) = (mn(n + 1) - S)(a^f + (m - lja^^ . 

(32) 

The instability intervals can be obtained imposing that 
the right-hand side of Eq. ([32]) be positive. The sign of 
C(„) is chosen so that the real part of Hk is positive if £ 
lies in an instability zone. 

If all the roots of the polynomial M( n \{z) are real and 
different, the definite integral in Eq. (|3l| ) can be ex- 
pressed in terms of elliptic integrals. Indeed, let z = l—y, 
and let the constants Pi, Di, i = 1, . . . , n, be defined by 



M (n) (l-y) 



1=1 



D; 



(33) 



The formula ( J3l| ) for the characteristic exponent p,k then 
reduces to 



^ * i—i 



(34) 



with 

nor 1 !*) 



n(/?rV) 



where IT(s|k) is the complete elliptic integral of the third 
kind || . 

We shall now compute the characteristic exponent for 
n = 1, 2, 3 using Eq. ([mJ). It shall be convenient to define 
the dimensionless momentum k = k/Mp. Using Eq. fl2C| ) 
and the definition of a, Eq. (Oh becomes 



£ = 



1 /2nD(m) ~ 9 , 

- ,\ ' k 2 + to + l)n(n + 1 

3 \ \a z 



(36) 



For n = 1, the Lame equation (^|) possesses two insta- 
bility zones, namely £ € (— oo, m) and £ € (1, 1 + to); 
see for instance Ref . |1| . The first instability zone is not 
admissible, since it would force the momentum fc to take 



imaginary values. Using Eq. (36), the second instability 
zone leads to the resonance condition 



\a 2 {l - 2m) 
2irD(m) 



<k 2 < 



\a 2 (m + 1) 
2tt£>(to) 



(37) 



Thus for to ^ i (i.e., for 5^0), there is an initial thresh- 
old for the resonant values of the momentum. If the x 
field had a mass m x , the instability interval (^) would be 
shifted by an amount —mt/Mp. For m x large enough, 
this could result in the disappearance of the threshold, 
or even of the whole interval. 

The polynomial Mm (z) is given by 



M (1) (z) = z + 



\-£ 



m 



Therefore 



Pi 1 = D 1 = T 



to — £ 



and 



r 2 

°(i) 



: (8 -m){8 -!){! + m- 8). 



Note that /3f 1 > 1 in the second instability zone. In 
Fig. |] we plot the characteristic exponent as a function 
of a and 10 12 fc 2 for a between and 10. The absolute 
maximum of the characteristic exponent for n = 1 is 
p,k — 0.147 at a — and k 2 = 1.08 • 10 -13 , in agreement 
with [pj . The maximum of the characteristic exponent for 
fixed a decreases monotonically with a. In contrast, the 
width of the resonance band of the squared momentum 
k 2 increases for a small, reaching its maximum at a = 
2.442. In order to compare the efficiency of the particle 
production with the pure A0 4 /4 model, we estimate the 
number density in units of Mp after 30 oscillations of the 
inflaton field (when the back reaction becomes significant 
in the A</> 4 model) as 



N(a) 



2tt 2 



~ k 2 e X2^ k K{n) d ^ 



(38) 



where / is the union of all the instability intervals. In 
Fig. ^| we plot the logarithm of the number density versus 



6 




.75 

10 12 fc 2 




N(a)/N(0) 



FIG. 6: Number densities ratio N(a)/N(0) for small values 



FIG. 4: Characteristic exponent as a function of a and W 12 k 2 
for n = 1. 



,JV(o) 




FIG. 5: Plot of log 10 iV(a) as a function of a for n = 1, 2, 3. 



<5. In Fig. [j] we represent the ratios N(a)/N(0) for small 
values of a. The number density for n = 1 (i.e., g 2 = A) 
is almost constant for a < 0.2. Therefore, in this case 
the effect of the term /z/(2</> 2 ) in the energy transfer from 
the inflaton field is negligible provided /x 1 / 6 < 10 _3 Mp. 
The number density decreases sharply for a > 0.5. 

We next examine the case n = 2. In this case the 
Lame equation (|^) possesses three instability zones, out 
of which the only relevant one is4 + m<£<2(l + m + 



\Jm 2 — m + l) , leading to 



3A<5 2 (2-m) 
2tt D(m) 



< k 2 < 



3\a 2 \/m 2 — m + 1 
irD(m) 



(39) 



The width of the resonance band of the squared momen- 
tum k 2 also increases in this case for small a, with its 

(2) 

maximum located at a — 0.858. The coefficients a\ 
(2) 

and a 2 defining the polynomial M(o\(z) are given by 



(omitting the superscript (2)) 

4-2to-£ (1 + m- £)(4 + m- £) 



ai 



3 m 



9m 2 



(40) 



One can immediately show that the roots of M( 2 ) (z) are 
real and different if £ lies in the above instability band. 
The coefficients /3i t 2 and D1.2 are respectively given by 



/3i,2 = -^(2 + 01 T \Ja\-Aa 2 ) , 



Dx.2 = 



±(2 + ai) + - 4 «2 



2(1 + ai + a 2 )\Ja\ - 4a 2 



(41) 
(42) 



Note that /3 1 2 > 1 in the insta- 



with ai : 2 given by 
bility zone. The coefficient C 2 2 ^ is given by 

1 



(1 + m - £ ) (4 + m - 8 ) (1 + Am - 8) 



(2) 81m 4 

x (£ 2 -4(l + m)£ + 12m) . 

The maximum of the characteristic exponent in the insta- 
bility band also decreases monotonically with a. The ab- 
solute maximum is = 0.036 at a = 0, k 2 = 7.67T0 -13 , 
four times smaller than the absolute maximum for n = 1. 
Correspondingly, the particle production is also much less 
efficient than in the case n = 1, see Fig. |^. The maximum 
particle production for n = 2 occurs when a = 0.245, see 
Fig. § 

The case n = 3 presents some unexpected effects. In 
this case, the Lame equation (|^) possesses four insta- 
bility zones, out of which the only relevant ones are 
£ e (4(1 + m),2 + 5m, + 2V4m 2 - m + l) and £ S 
(5 + 2m + 2\/m 2 - m + 4, 5(m + 1) + 2V4m 2 - 7m + 4) . 
In terms of the dimensionless momentum k, these are 



< k 2 < 



3Ac 



27rD(m) 



(m - 2 + 2 V4m 2 - m + l) , (43) 



7 



and 



3A 



2nD(m 



a 2 i 

—— (l - 2m + 2y/m 2 - m + 4) < fc 2 



< 



_3Aa 2 _ 
2irD(m) 



(l + m+ 2\J\m 1 - 7m + 4) . (44) 



For both resonance bands, the width increases for small 
<S, with the respective maxima located at a = 1.052 and 
a = 0.533. The coefficients of the polynomial M( 3 )(z) 
defined in Eq. (ROT) are given by 



ai 



9 — 6m — £ 
5m 



2£ 2 + (4m - 26)£ + 27m 2 - 51m + 72 
« 2 = " ^2 ( 4 5) 



75m 



a3 = 



(f 3 - 2(4m + 7)£ 2 + (16m(m + 5) + 49)£ 

5TO J V 



225m 
- 12(m+ l)(8m + 3) 



The coefficients C? 3 n and Di, i = 1,2,3, are then 
easily obtained from Eqs. ( |32"1 ) and (|33|). The resulting 
expressions are very cumbersome and shall not be dis- 
played here. It may be shown that the coefficients j3i are 
all real and different in the resonance bands. 

Just as in the previous cases n = 1,2, the maximum of 
the characteristic exponent decreases monotonically with 
a for both resonance bands. For the lower resonance 
band given in Eq. ([43|), the absolute maximum is fj,k = 
0.159 at a = 0, fc 2 = 2.09- 10~ 13 , while for the higher one 
in Eq. the maximum value is fik — 0.0078 at a = 0, 
k 2 = 1.857 • 10~ 12 . The particle production at a = is 



two orders of magnitude more efficient than in the case 
n = 1. However, it decreases with a much faster than in 
the previous cases n = 1, 2, see Figs. 0, ^. 



IV. CONCLUSIONS 

In this paper we have characterized the most general 
scalar potential for the inflaton field leading to the Lame 
equation for the matter field modes in a Minkowskian 
background. The resulting potential possesses a term of 
the form (j,/(2(j) 2 ) in addition to the terms A^ 4 /4 + if (j> 2 /2 
already studied in the literature. We have analyzed the 
effect of this new term in the preheating era after infla- 
tion in the particular case K = 0. Exact expressions 
for the resonance bands and the characteristic exponents 
have been derived for certain values of the coupling con- 
stant between the inflaton and the matter fields. The 
effect of the new term in the particle production is vir- 
tually negligible provided z^ 1 / 6 < 10~ 3 Mp, even though 
the inflaton potential is modified in an essential way near 
the origin. However, for /i 1 / 6 ^> 10~ 3 Afp matter produc- 
tion is heavily suppressed by the new term as compared 
to the pure A</> 4 /4 model. The situation in this respect 
is expected to remain unchanged for other values of the 
coupling constant. 
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